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ABSTRACT 


A mathematical model of the mass characteristics of composite bodies has been 
programmed for digital computer processing. The description and equation development 
of the mathematical model are presented. The technique is of particular interest for 
spacecraft whose size or mass prohibits the use of experimental measurement in an in- 
ertia laboratory. The computer program is written in the FORTRAN IV (version 13) 
language and is operational on the IBM 7094 Direct Couple System of the Lewis Research 
Center. A sample problem illustrates the input data requirements and the output format 
of the program. The use options of the program and program limitations are fully de- 
scribed. 



A FORTRAN CODE FOR COMPUTING THE PRINCIPLE MASS 
MOMENTS OF INERTIA OF COMPOSITE BODIES 
by James E. Cake 
Lewis Research Center 

SUMMARY 

A mathematical model of the mass characteristics of composite bodies has been pro- 
grammed for digital computer processing. The description and equation development of 
the mathematical model are presented. The technique is of particular interest for space- 
craft whose size or mass prohibits the use of experimental measurement in an inertia 
laboratory. 

The computer program is written in the FORTRAN IV (version 13) language and is 
operational on the IBM 7094 11/7044 Direct Couple System of the Lewis Research Center. 
A sample problem illustrates the input data requirements and the output format of the 
program. The use options of the program and program limitations are fully described. 

INTRODUCTION 

Studies of the rotational dynamics of spacecraft require the knowledge of the mass 
characteristics of the spacecraft. These characteristics include spacecraft weight, 
location of the center of mass, principal mass moments of inertia, and the directions of 
the principal axes of inertia. Laboratories to experimentally measure these character- 
istics may not be readily available or cannot handle spacecraft of large size and mass. 
Therefore, a mathematical model which treats the spacecraft as a composite body is 
desirable. 

This report describes the mathematical model which has been developed and pro- 
grammed for digital computer processing. The computer program is written in the 
FORTRAN IV (version 13) language and is operational on the IBM 7094 11/7044 Direct 
Couple System of the Lewis Research Center. 

The report places particular emphasis on data preparation for the computer pro- 
gram. A sample problem is included for both illustrative purposes and as a verification 



of the mathematical model. The program has been used to calculate the mass character- 
istics of the Space Electric Rocket Test vehicle (SERT II). 


THE PROBLEM 


The mathematical model of the mass characteristics of a composite body is based on 
the fact that the moments and products of inertia of composite bodies are found by sum- 
ming the moments and products of inertia of the individual components or particles. This 
may be shown from a few fundamental principles of rigid body and particle dynamics. 

From reference 1, the components of the angular momentum of a rigid body about a 
fixed point of rotation, 0, are: 


\ =Ix 
x 0 



CO - P CO - P CO 
0 x 0 xy 0 y 0 xz 0 z 0 

CO + I CO - P CO 
yx 0 x 0 y 0 y 0 yz 0 z 0 

CO - P CO + I CO 
zx 0 x 0 zy 0 y 0 z 0 z 0 






(1) 


where 1,1. and I are the mass moments of inertia and P„_ , P„ Y , etc. , are 
^0 y 0 z 0 y 0 z 0 
the mass products of inertia of the body. 

The total angular momentum of a system of particles is the sum of the angular mo- 
menta of the individual particles. 
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Substitute equation (1) into equation (2). By expanding the right side, the compo- 
nents of the angular momentum of the system can be expressed in a form similar to 
equation (1). The moments and products of inertia of the system of particles or compos- 
ite are found to be the sum of the moments and products of inertia of the individual com- 
ponents . 

Description of Mathematical Model 

The spacecraft or other rigid body is composed of n discrete components of mass 
n^, each with its center of mass at (x^ y^ in the reference axis system defined by 
figure 1. 

Definition of the center of mass . - The center of mass (x, y, z) of the body with re- 
spect to the reference axis system is defined as 

i>i*i 



n 

i=l 

E-v. 

i=l 

i-ft 



n 

Z m i 

i=l 

Definition of the inertia tensor . - The mass moments of inertia of each component 

are L. , I„ , I which are defined about axes parallel to the reference axis system. 

*i y i z i 

The center of mass axis system as defined in figure 1 has its origin at the system center 
of mass and is parallel to the reference axis system. Then the system mass moments of 
inertia about the center of mass axis system are: 
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n 

X is + m ‘t <j ■ z ‘ >2 + (y ■ y ‘ )2 ]} 

i=l 

X { Iy i + mi [ <Z ' Z ‘ )2 + <Z " X ‘ )2 ]} 

i=l 

n 

X pi + “i [ (S - */ + ( y ■ y i )2 ]} 


The mass products of inertia defined about axes parallel to the reference axis sys- 
tem are P P , P P , and P P . The system mass products of inertia 
xy i yx i yz i zy i xz i zx i 

parallel to the reference system are 


P — = p — = 

xy yx 


n 

Ih (x - X i )( y - y t ) + P xy J 


P — = p — = 
yz zy 


n 

^(y - Yi)(z - Z t ) + Py Z . 


P = P = 

XZ Z X 


Y_, h <x ■ x i )<z ■ z ‘’ + % 
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The inertia tensor is defined as 


_I 11 

h2 

*13 


' *x 

- P xy 

- P xz 

*21 

X 22 

*23 

= 

“ p yx 

■y 

1 

‘Cl 

N| 

hi 

X 32 

*33 


- p zx 

■ P zy 

% 


From reference 2, the principal mass moments of inertia are found from the determinant 


11 - 1 ) 

*12 

*13 

*21 

d 2 2 - *> 

*23 

*31 

*32 

(*33 " 


Evaluating the determinant yields a cubic equation in I whose roots are Ip ^ I3, the 
three principal mass moments of inertia. 

Orientation of the principal axes . - If the products of inertia are nonzero in the in- 
ertia tensor, then as shown in reference 2, the following method can be used to solve for 
the directions of the principal axes. 


U) co 

(I 22 - i)_Z + i M _* = -i 


23 

(O U) 

X X 


21 


(?) 


0 ) 


co. 


>32 f + (I 33 - *> ^ - - 1 

X 5 


31 


( 8 ) 


Substitute each principal inertia I- for I in equations (7) and (8). Solve for co /w 

J y j/ X j 

and co . By arbitrarily choosing o> = 1, then co„ /o> = w r and c o„ /co = co . 

V j j y j/ x j y j . V X j Z j 

Then the magnitude of the angular velocity vector 


r. = */w 2 + co 2 + co 2 = 

] M x i y i z i V 


. 2 2 

1 +( V + C V 
J 
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The direction cosines are found to be 


* A x = — 
A l x r x 


t A ^ t. A v 


A 2 x r 


A 3 x r. 


U) U> 0) 

l -Jl l --ll l — J3 
V " r, V ' r n V r 0 


co„ 


^ A l z r< 


u> 


1a 2 z r r 


u> 


^ A 3 z r. 


where Aj is the axis about which Ij is defined, etc. 

The set of direction cosines represents the rotation matrix between the center of 
mass axis system and the principal axes of the body. 


V 


l A 1 x 

l A x y 

^AjZ 


C 

X 

< A 2 

> = 

1a 2 x 

1a 2 y 

' A 2 z 

< 

y l 

A 3 
C J 


! A 3* 

* a 3* 

' A 3 z . 


z 

^ J 


(9) 


If four products of inertia in the inertia tensor are zero, an alternate method is re- 
quired to solve for the direction cosines. Consider the case where x is a principal 
axis. Then the inertia tensor is 



x n 

0 

0 



0 

*22 

*23 

(10) 


0 

*32 

CO 
CO 
t— 1 
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The direction cosine matrix between the center of mass axis system and the principal 
axes represents a single rotation a in the yz plane. 


[l] = 


0 cos a sin a 


[_0 -sin a cos oj 


(ID 


From reference 2, 


[i'l = M [i] [J] T 


( 12 ) 


where [I* ] is the inertia tensor about the principal axes of inertia. Substituting (8) and 
(9) into (10) and solving for I^, 


2 1 2 
Ig 2 = I32 cos a + - (I33 ■ I 2 2) sin 2a ■ I 2 3 sin a 
2 


and since I 32 = ^3, then 


I32 = I32 cos 2a + - (I 33 - I 2 2) sin 2a 
2 


In the principal inertia tensor I 32 = 0. Therefore, reducing equation (11) yields 



2I 32 

T 22 " l 23 


A similar approach is used to find the direction cosines whenever y or z is a principal 
axis. 


Accuracy of Mathematical Model 

The accuracy of the mathematical model is dependent on the accuracy of the mass 
characteristics of the individual components. Ideally, accurate mass characteristics of 
each component would be found in a laboratory which measures the quantities experimen- 
tally. If this is not possible, the component may be assumed to be homogeneous in mass 
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distribution and of a regular geometric shape. Standard formulae may be used to calcu- 
late the moments and products of inertia about axes parallel to the reference system. If 
the center of mass is not found experimentally, assuming a homogeneous mass distribu- 
tion locates the center of mass at the geometric center of each component. 

Less accurate results will be obtained if any of the above assumptions or simplifica- 
tions are made in calculating the inertia tensor and center of mass rather than experi- 
mentally measuring these mass characteristics. For any method, the resulting compo- 
nent inertia tensor would be the input data to the computer program. 

DESCRIPTION OF DIGITAL COMPUTER PROGRAM 

The computer program is written in the FORTRAN IV (version 13) language and is 
operational on the IBM 7094 11/7044 Direct Couple System of the Lewis Research Center. 

The main program, READ, controls the input and output functions desired by the 
user. Subroutine INRTIA calculates the principal mass moments of inertia and the direc- 
tions of the principal axes. Subroutine OUT lists the input data for the individual compo- 
nents, the system inertia tensor, the principal inertias and directions of the principal 
axes, the system weight, and the system center of mass. Schematic flow diagrams of 
READ, INRTIA, and OUT are presented in figures 2 to 4, respectively. Subroutine 
EIGEN is called by INRTIA to calculate the eigenvalues (principal inertias) of the system 
inertia tensor, and subroutine CUBIC solves a cubic equation. 

Common statement symbols used in the program are defined in appendix B, and a 
complete program listing is given in appendix C. An understanding of the input data 
preparation and output data interpretation is essential to using the program. 


Program Use Options 

The input and output portions of the program have been written to make its use sim- 
ple and efficient. The initial input data to the program is placed on punch cards. This 
data may then be written on a magnetic tape. For future calculations, the data is input 
using the tape with any changes to the data array made by card input. The tape may then 
be changed to the new data array or retained as the original data set-up. 

The input options, which are assigned the FORTRAN variable name NOPT, are: 

(1) NOPT = 1; the data array is input from cards. 

(2) NOPT = 2; the data array is input from magnetic tape. No changes or additions 
are made to the array. 

(3) NOPT = 3; the data array is input from magnetic tape. Changes or additions are 
made to the array . 

(4) NOPT = 4 ; the data array is not input through READ. 
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The output options, assigned the variable name NOUT, are: 

(1) NOUT = 1; subroutine OUT is used for output data. 

(2) NOUT = 2; subroutine OUT is not used. 

The magnetic tape options, assigned the variable name NTAPE, are: 

(1) NTAPE = 1; the new data array will be written on tape. 

(2) NTAPE = 2; the tape is left unchanged from the previous data array. 

Program Input Preparation 

There are eight possible types of input cards to the program. The use of the card 
will vary with the program use option. 


Card 

FORTRAN 

Format 

Card 

Description 

Units 

type 

1 

variable 

NOPT 

15 

columns 

1-5 

Input option 



NOUT 

15 

6-10 

Output option 



NTAPE 

15 

11-15 

Tape option 


2 

DESCRI 

5A6 

1-30 

Component description 



NCOMP 

15 

31-35 

Component number 



IP ART 

15 

41-45 

Part number 



MODEL 

15 

51-55 

Model number 



IDATE 

16 

60-61 

Month 


3 

NCOMP 

15 

62-63 

64-65 

1-5 

Day 

Year 

Component number 



MASS 

E12.4 

6-17 

Component weight 

lb 


IXC 

E10.3 

18-27 

Vi 

slug-ft z 


IYC 

E10.3 

28-37 

slug-ft 2 


IZC 

E10.3 

38-47 

S 

x i 

slug -ft 2 


XC 

E10.3 

48-57 

inches 


YC 

E10.3 

58-67 

yi 

inches 


ZC 

E10.3 

68-77 

z i 

inches 

4 

PXYO 

E10.3 

1-10 

% 

slug-ft 2 


PXZO 

E10.3 

11-20 

P xz 

slug-ft 2 


PYZO 

E10.3 

21-30 

V, 

slug-ft 2 
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Card 

FORTRAN 

Format Card 

type 


columns 

5 

IADD 

11 1 


N 

14 2-5 

6 

Same format as card 2, except all 

7 

Same format as card 3, except all 

8 

Same format as card 4, except all 


Description Units 

Change or addition indicator 

List number 

zeros punched in the entry fields 

zeros punched in the entry fields 

zeros punched in the entry fields 


Each component requires a set of three input cards (cards 2, 3, 4). For any compo- 
nent input by cards, a set of "zero field" data (cards 6, 7, 8) must be the last set read 
by the program. 


Data Array Change Using Tape Input 

The program assigns a list number to every component in the data array including 
the "zero field" data. When changing the data array after a tape input the following 
quantities must be specified on card 5: 

(1) IADD = 0; the data of a component in the array is being changed. 

= 1; a new component is being added and the length of the array is increased. 

(2) N; the list number of the data to be changed or added. 

When adding components, the list number of the first addition is assigned the list number 
of the "zero field" data which was on tape. For either component additions or changes, 
the new set of "zero field" data required to end the read phases is always considered an 
addition to the list. 

The following examples illustrate the card set-up for the various input options: 

Option 1: Cards 1, 2, 3, 4, 2, 3, 4, 2, 3, 4, 6, 7, 8 
Option 2: Card 1 

Option 3: Cards 1, 5, 2, 3, 4, 5, 6, 7, 8 
Option 4: Card 1 

Batches of inertia calculations may be run using the program. 

Appendix D contains a sample problem illustrating the input data requirements and 
the output format of the program. 


Program Limitations 

Because of core limitations, the program has reserved storage locations for a com- 
posite body with a maximum of 850 components. For larger composites the inertia cal- 
culation should be made as separate jobs. 
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The magnetic tape is programmed to operate from logic unit number eight. There- 
fore, any control cards preceding the program must specify the tape to be "set-up” on 
that unit. 


CONCLUDING REMARKS 

The purpose of this report has been to describe a computer program capable of com- 
puting the mass characteristics of composite bodies. The input words and output data 
have been fully defined by the description of the mathematical model. A sample problem 
was used to illustrate the input options, input data preparation, and output data interpre- 
tation.' 

The accuracy of the mathematical model may be increased by experimental measure- 
ment of the mass characteristics of the component. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, December 6, 1968, 

704-01-00-17-22. 
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APPENDIX A 


MATHEMATICAL SYMBOLS 


A l’ A2, A 3 

H x ’ H v ’ H z 
x 0 y 0 z 0 




I 

l 

m 

P 

r 


x i’ y P z i 


x,y, z 


a 

U> 

Subscripts: 

x 0 ’ y 0’ z 0 

x,y,z 

x',y',z’ 

Superscript: 

T 


coordinate axes defining the principal mass moments of inertia 
angular momentum of a system of components about x Q , y 0 , and 

2 2 \j yj 

respectively, slug -ft /sec; kg-m /sec 

angular momentum of a component about Xq, y^, and Zq, respectively, 
slug-ft^/sec; kg-m^/sec 

2 2 

mass moment of inertia, slug -ft ; kg-m 
direction cosine 
mass, slug; kg 

2 2 

mass product of inertia, slug-ft ; kg-m 

magnitude of the angular velocity vector, rad/sec 

component center of mass in the reference axis system, ft; m 

coordinates of the system center of mass in the reference axis system, 
ft; m 

angle of rotation, deg 
angular velocity, rad/sec 

reference axis coordinate system 

center of mass coordinate system 
principal axes coordinate system 

matrix transpose 
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APPENDIX B 


COMMON STATEMENT SYMBOLS 

A element of inertia tensor 

DCSN(J, J) direction cosines 

DESCRI(J, J) description of component 

IDATE(J) date of inclusion into inertia calculation 

IPART(J) part number of component 

KC(J) moment of inertia about the x axis of the component 

IYC(J) moment of inertia about the y axis of the component 

IZC(J) moment of inertia about the z axis of the component 

IXB(J) component moment of inertia translated to the x axis 

IYB(J) component moment of inertia translated to the y axis 

IZB(J) component moment of inertia translated to the z axis 

KY(J) component product of inertia translated to the xy plane 

IXZ(J) component product of inertia translated to the x z plane 

IYZ(J) component product of inertia translated to the yz plane 

MASS(J) component weight 

MODEL(J) component model number 

NCOMP(J) component number 

PI1, PI2, PI3 principal mass moments of inertia 
PXYO(J) product of inertia about the xy plane of the component 

PXZO(J) product of inertia about the xz plane of the component 

PYZO(J) product of inertia about the yz plane of the component 

SWEIGH total system weight 

XBAR system center of mass distance from Xq axis 

YBAR system center of mass distance from yp axis 

ZBAR system center of mass distance from Zq axis 

XC(J) component center of mass distance from Xq axis 
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YC(J) 

ZC(J) 

XWRTCG(J) 

YWRTCG(J) 

ZWRTCG(J) 


component center of mass distance from Yq axis 
component center of mass distance from Zq axis 
distance from component to x axis 
distance from component to y axis 
distance from component to z axis 
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APPENDIX C 


FORTRAN LISTING 


$ I B FTC R F AO 


C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 


c 

c 

r. 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


USF OPTIONS 

NOP T = l READ INPUT FROM, CARDS 

NOPT = 2 READ INPUT FROM TAPE, NO CHANGES 

NOP T = 3 READ INPUT FROM TAPE. MAKE CHANGES AND ADDITIONS 

NOP T = 4 INPUT THROUGH SUBROUTINE READ IS NOT USED 

NOti T = 1 USF OUTPUT SUBROUTINE OUT 

NQUT = 2 SUBROUTINE OUT IS NUT USED 

N T A P H = 1 NEW DATA ARRAY WRITTEN ON TAPE 

NTAPE=2 TAPE LEET UNCHANGED FROM PREVIOUS CATA ARRAY 


R EM ARK. S 

IF THE INPUT AMD OUTPUT ROUTINES ARE NOT USED, CATA TO 
SUBROUTINE INRTIA IS TRANSMITTED THROUGH 

COMMON /INPUT/ 

COMMON /OUT PTA 

COMMON /TENSOR/ 


DESCRIPTION OF AXES 


1. TEE REFERENCE AXIS SYSTEM IS DEFINED BY THE USER 

2. THE CENTER OF MASS AXIS SYSTEM IS PARALLEL TO THE REFERENCE 
SYSTEM WITH ORIGIN AT THE SYSTEM CM 

3. INPUT MOMENTS AND PRODUCTS OF INERTIA ARE ABOUT AXES 
PARALLEL TCI THE REFERENCE AXIS SYSTEM WITH ORIGIN AT THE 
COMPONENT C.M. 


DESCRIPTION OF PARAMETERS 


MAS S 
XC.YC.ZC 


WEIGHT OF COMPONENT (LBS) 

LOCATION OF COMPONENT IN REFERENCE 


AXIS SYSTEM (INCHES) 

IXC, IYC, IZO. MOMENTS OF INERTIA CF CCMFONENT IN 

COMPONENT AXIS SYSTEM (SLUG-ET2) 
PXYO.PXZU.PY/.O PRODUCTS OF INERTIA OF CCMFCNENT IN 

COMPONENT AXIS SYSTEM ISLUG-FT2) 
XBAR.YBAR.Z8AR CENTER (JF MASS (INCHFS) 

XWR TOO. YWRTC.G, ZWRTCG DISTANCE OF COMPONENT FROM C.M. (INCHES) 
'IXrt.IYB.I7Ii MOMENTS OF INERTIA CF COMPONENT IN C.M. 


ixy, rxz, IYZ. 

SwF IGH 

P ft .PI2.PT 3 

COST) 

A 


AXIS SYSTEM ( SLUG-FT2 ) 

PRODUCTS OF INERTIA OF COMPONENT IN C.M. 
A X I S SYSTEM ( SLUG— FT2 ) 

WEIGHT OF SYSTEM (LBS) 

PRINCIPAL MOMENTS OF INERTIA OF SYSTEM 

DIRECTI ON CCSI NE SET 

ELEMENT OF SYSTEM INERTIA TENSOR 
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COMMON/ INPUT/ OESCRI ( 850,5) ,NC OMPI 850) , I PART {850) , MODEL (850 ) , 

♦ I CATE! 850) .MASS! 8 50 ) . I XC I 850) . I YC I 850) ,1 ZC (850 > ,XC 1850 ) ,YC I 850 ), 

♦ ZC ( 85 Ol.PXYOl P 50) .PXZOI850) ,PYZ0(850) 

COMMON /OUTPT/XWRTCG I 85 0 , YWR TO 0(850) , Z WRTCG < 85 0 ) , I XB ( 3 50 ) . I YB< 850) 

♦ . IZBI 8 50 ) . IXYI 3 50) , I XZ ( 85 0) , I YZ (850) .COSO (3 ,3) ,P I 1 , FI 2, PI 3, SWE IGH, 
♦XFAR.YBAR.Z8AR 

COMMON/TENSOR/ A 1 1.A12.A13 .A22.A23.A33 
ECU TVAI. ENCEI A 12.A2 1 ) . ( A 13, A 31 ) , (A23, A3 2 ) 

REAL MASS, IXC, I YC, I ZC 

4 I C N T = 0 
REWIND .3 

C SELECT INPUT AND OUTPUT OPTION 

READ(5«103) NOP T.NOUT .NTAPE 
GO TO ( 1, 2. 2, 500) ,NUPT 
C CARD INPUT 

1 ICNT=TCNT+L 
N = I CNF 

K EAD15. 101 KITE SCR I ( N , J ) , J =1 ,5 ) , NC GMP < N ) , I PART (N) , MODEL ( N) , I DAT E( N) 
R E ADI 5, 10 0) NCOMP ( N ) , MASS ( N) »I XC(N) .1 YCIN) ,1 ZC(N) ,XC(N),YC(N),ZC( 
*N ) 

READ! 5, 1C2 ) PXYO IN ) , P X ZU I N ) .PYZQIN) 

I F < MASSIN) .NE .0.0) GO TO 1 
GO TO 500 
C TAPE INPUT 

2 ICNT=lCNT+l 
N = ICN r 

READ! 8)1 OESCRI (N,J) , J = I , 5 ) .NCGMPIN) , I PARTI N) .MODEL IN), I DAT El N) , 

♦ MAS SIN). IXC(N), I YCIN), IZC(N).XCIN) .YCIN) ,ZC IN) , PXY C I N) , PXZ 01 N.) , 
♦PYZOIN ) 

IF! MASS! N ) .NE .0. 0) GO TO 2 
I FI NOP T .EQ.2 ) GO TO 500 
I CN 1= l CNT-1 
10 R FAD! 5, 104 ) IADD.N 

READ! 5, 10 1 ) IDE SCR I (N,J),J = 1,5> .NCOMP IN) , I PART I N) .MODEL IN) , I CAT E{ N) 
R E ADI 5, 100) NCOMP! N > .MASS ( N ) , I XC ( N ) , I YC I N) , I ZC I N ) , XC ( N ) , YC I N ) , 

* Z Cl N ) 

R EADI 5, 102 I PXYOIN ) ,PXZO IN) .PYZOIN) 

I FI I ADD . EQ . 1 ) ICN T=ICNT+1 
I FI MASS! N) .NE .0.0) GO TO 10 
500 CALI INKT I A ( ICNT) 

C SELECT OUTPUT OPTIONS 

I FINTAPE .EQ .0) GO TO 6G0 
REWIND 8 
DO 5 N =1 , I CNT 

5 WRI T El 8) I QESCR I IN. J ) . J = 1 ,5) .NCOMP! N) , I PART IN) .MODEL IN) , l DATE! N ) , 

♦ MASS IN ). IXC IN ) ,1 YCIN) . I ZC I N ) ,XC(N) ,YC IN) .ZCI M , PXYOIN) , PXZOI N) , 

♦ PYZOIN ) 

600 IFINOlJT.EQ.l ) CALL OUT! ICNT) 

GO TO 4 

100 FORMAT! 15, E 12.4, 6F 10.3 ) 

1.0 I FORMAT I 5A6, 15. 5X. 1.5 , 5X.I 5.4X.I 6) 

102 FORMAT! 3E 1C. 3) 

103 FORMAT I 15, 15.15) 

104 FORMAT I 11,14) 

FND 
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1IRFTC INERTIA 


SL8P0UTINF INRTIA(ICNT) 

C DM MON / I NPli I / DESCRK F.50,5) .NC0MPIH50) ,1 PART (850) , NQDELI850), 

* I CATE( 8 50 ) .MASS( 8 50) , I XC ( 85 0) , I YC ( 85 0) , I ZC (850 ) ,XC(850) ,YC(850 ) . 

*= 7C( B50) ,PXYO( 850) .PXZG(85C) ,PYZfi(850) 

CQMMOM/OL TPT/XWRTCG ( 85 C> . YKR TC G ( 85 0 ) , ZWRTCG(350) , IXB (a 50 ) , I YB( 850) 
!7_R( 850) . I XY { 850) . I XZ ( 850) .1 YZ (350) ,CCSD(3 ,3) ,P II , F 1 2 , PI3.SWE IGH, 
*XFAR,YBAR.ZRAR 

GO M MON /TENSOR/ A I 1 . A12, A13 ,A22, A23, A3 3 
D JMENS ION PI(3).Z(3)*Y(3) 

EQUIVALENCE ( PXY.P YX) , ( PYZ, PZY ) , ( PXZ, PZX) 

ECU IVAL FNCE ( A 1 2, A 2 1 ) , ( A 23 , A32 ) . ( A 13 , A 31 ) 

ECU I VALENCE (All. ROLL ).(A?2»PI TCH) , (A3 3, YAW) 

ECU IVAL FNCE ( P I ( l ) ,P I 1) 

R F AL MASS . IXC. I YC , I ZC , I XB ♦ I YB . I ZB , I X Y , I XZ , I YZ 

r. 

C LOOP FOR CALCULATING SYSTEM CENTER OF MASS 

C 

SXCM = 0 .0 

S YCM=0 .0 

S 7 C M = . ) . 0 

S M A S S = 0 . 0 

DC 5 N = 1 . ICNT 

M A S S ( N ) =MA S S ( N ) /3-2 . 1 74 C-'« S 

X C ( N ) = XC ( M ) / 1 2 . 0 

YC(N)=YC(W)/12.0 

Z C ( N ) = ZC(N ) Z12 .0 

X f M = XC ( N )-MASS(N) 

YCM=YC( N )*MASS (N) 

Z CM = ZC ( N )*MASS(N) 

SMASS=SMASS+MASS(N ) 

SXCM=SXCM+XCM 
SYCM=SYC M + YCM 
5 S7CM =SZCM+ZCM 
X R A R = .S XC M / SMA S S 

Y FAR-SYCM/SMASS 
7 RAR = S7CM/SMASS 
N = ION T 

c 

r LOOP FOR CALCULATING SYSTEM MASS MOMENTS AND MASS PRODUCTS 

C 

R EL L=U .0 

Y AW = 0 .0 

P I T C H = 0 » 0 
P X Y = Q . 0 
P XZ =0.0 
P Y 7 = 0 . 0 
00 3 K = 1 ,N 
X = X fi AK- XC ( K ) 

Y =Y 8 AR-YCI K ) 

7 =Z BAR— 7 C ( K ) 
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n o n n n n n o n n non 


I X ft I K ) = I XC ( K ) +M A SS IK)* (Z*Z+Y*Y) 

IYB(K) =1 YCIK 1+MASS (K)*( Z*Z+X*X ) 

I Z B ( K) =IZC(K 1+MASSIK1* ( X$X+Y*Y) 

R CL L =ROL L + I X8 I K ) 

P TTCH=P ITCH+IYBIK) 

Y AW=Y4W+ IZ8 (K ) 

TXY (K >=MASS(K )*X*Y+PXYU(K) 

IXZ (K> =MASS(K )*X*Z+PXZO(K) 

IY7(K)=MA$S(K)*Y*Z+PYZ0(K) 

P XY=PX Y+ I XY ( K ) 

P YZ = PYZ+ IYZ (K ) 

PXZ=PXZ+ I X 7 { K ) 

S FT UP INERTIA TENSOR. OTHER ELEMENTS BY EQUIVALENCE 

A 1 2=-P XY 
A23=-PYZ 
A n=— PX7 
3 CONTINUE 

I F ( A12.E0.0.0.AND.A13.EQ.0.0.AN0.A23.EQ.0.0) GC TO 24 
IF< A12.E0.0.0. AND. A13.EQ. 0.01 GO TO 21 
IF(A12.EQ.0.0.AND.A23.EQ.0.C) GO TO 22 
IF! A13.EQ .O.O.AND.A23.FQ.O.C) GO TO 23 

SOLVE THE EIGENVALUE PROBLEM FOR THE PRINCIPAL INERTIAS 

CALL E IGEN( PI1,PI2,PI3) 

I F ( P Il.EO.O.O) GO TO f 

DERIVE DIRECTION COSINE MATRIX 

DO <5 1 = 1,3 

Z (I > = ( — A 2 t*A32+A3L*< PITCH- PI (I) ) > / (A23*A32— I PI TCH-PI ( I ) ) * 

II YAW-P I I III) 

Y ( I > = ('-A21-A23*ZI I )) / (PI TCH-PI ( I ) ) 

R=SORT( 1 .O+YII )*Y{ I ) + zm*zm ) 
cosni i,n = i .n/R 

COSCI I,2)=Y( I)/R 
COSO! 1 , 3 ) = zm /R 
9 CONTINUE 
GO TO 91 

DIRECTION COSINES AND PRINCIPAL MOMENTS WHEN XX IS A PRINCIPAL 
AXIS. 

21 ALP HA= . 5* ATAN( I2.0*A23) /< A22-A33) ) 

S A = S I N ( ALPHA ) 

CA=COS{ ALPHA ) 

COSDI l , 1 ) = 1 .0 
COSCI 2, 1 ) = 0 .0 
CCSDI 3,11=0.0 
CCSDI l , 2 ) =0.0 
COSCI 2, 2 1=CA 
COSO (3,2 ) =— SA 
COSDI 1,3 1 = 0.0 
CPSC!2.3)=SA 
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ccsni 3,3 )=ca 

PI1=AI1 

P I2=A2?*CA*CA+A33*SA*SA*A23*SI N(ALPHA*2.0) 

P I 3=A22*SA*SA *A3 3- CA*C A- A 2 3* SI N( ALPHA* 2. 01 
f.C TO 91 
C 

C DIRECTION COSINES AND PRINCIPAL MOMENTS WHEN VY IS A PRINCIPAL 

C AXIS 
C 

22 ALPHA-. 5* AT AN I I 2 . 0* A 1 3 ) / < A3 3- A1 1 )) 

C A=COS( ALPHA ) 

SA = S IN ( ALPHA ) 

ccsni i, n=cA 
COSOI 2, I >=0.0 
COSOI 3. 1 ) = SA 

ccsni i.« 2 ) = 0 .o 

Cf SDI 2. 2 > = 1 .0 
ccsni 3. 2 ) = 0 . 0 
ccsni 1 , 3 ) =— SA 

C0SDI2. 3 ) = 0 .0 
COSOI 3. 3 ) =C A 

P I 1=AL l*CA*CA*-A3 3* SA* SA- A 1 3*SI N(2. C* ALPHA) 

P I 2 = A22 

P I 3 = A 1 1*SA*SA*A33*CA*CA + A13*SI M2. 0* ALPHA) 

CO Tn 91 

c 

C n IP. EC r ION COSINES ANO PRINCIPAL MOMENTS WHEN ZZ IS A PRINCIPAL 

C AXIS. 

C 

23 ALP HA= . 5* AT AN I I 2 . 0* A 1 2 ) / I A1 1-A22 ) ) 

C A = COS I ALPHA ) 

SA=S INI ALPHA ) 

COSOI 1< 1 )=CA 
COSOI 2 . 1 )=-SA 

ccsni 3. l ) = o.o 
ccsni l, 2 )=SA 
ccsni 2 . 2 ) = C A 
C0S0I3.2 ) = 0.0 
COSOI l. 3 ) = 0.0 
COSOI 2. 3) =0.0 
C C S n I 3 , 3 ) = 1 . 0 

P I I = A l 1*CA*CA + A22*SA*SA*A12*SIN(2. 0*AL PH A) 

P 1?=A1 1*SA*SA+A22*CA*CA— A12*SI M2. 0* ALPHA) 

P I 3 = A 3 3 
GO TO 9 1 
C. 

C DIRECTION COSINES AND PRINCIPAL MOMENTS WHEN THE PRINCIPAL AXES 

C ARE ALIGNED WITH THE CENTER OF MASS SYSTEM 

c 

2A cnsni i , n= i -o 

COSOI 2. 1 >=0.0 
COSOI i. 1 )=0.0 

ccsni i. 2 i = o. c 

C.CSOI 2.21 = 1.0 
'COSOI 3. 2 1 = 0.0 
f.C SOI 1 . 3 ) = 0 .0 
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CGSD ( 2. ? )=0.0 
CDSUl 3.3 )=1 .0 
P 1 1=A 1 l 
P I2=A 22 
PI3 = A33 
GO TO 91 
C 

C CCNVER SIGNS FOR OUTPUT BLOCK 

C 

91 5>V»FIGrl-=SMASS*32. 1 74 04 9 
Xf AR=XBAR*12.0 

Y EAR =Y EAR* 12 .0 
7 P AR = 7 R AR* 1 2 . 0 
OC 6 K=1,N 

MASS1K. )=MASS(K )*32. 174049 
X C ( K ) = XC ( K )* 1 2 . 0 
YC(K ) = Y 0. ( K )* 12 . C 
7 C 1 K ) = ZC ( K ) * l 2 . 0 
XWRTCG! K )=XRAR— XC( K ) 

Y WR TCG< K ) = Y8AR- YC ( K ) 

6 ZVsRTCGIK )=ZBAR-ZC(K) 

GH TO A 
8 WRITE! 6. 233 ) 

23.3 FORMAT ( 9 3H0EX ECUTI ON HAS BEEN TERMINATED. SOLUTION TO THE EIGENVA 
ILliE PROBLEM RESULTED IN COMPLEX ROOTS) 

A CONTINUE 
RETURN 
END 
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tIBFTC PUT 


SI.RRDU1 IMF OIJTIICNT) 

COMMON/ INPUT/ RE SC R I ( 8 50 . 5 ) , NC GMPC 850 ) , I P AR T ( 8 50 ) , MODEL (850 ) , 
l CAT E( 8 50 > .MASS! 6 50 ) , 1 XC I 65 0) , I VC ( 850) , l ZC ( 850 ) , XC ( 8 50 ) , YC ( 8 50 ) , 

* ZC! 850 ) , PXYfH 8601 .PXZD ( 85 GJ , P YZO( 850 ) 

CCMMUN/OiUTP T /XVJR TCG ( 850) , YWRTCG ( 850) , ZWRTCG (850) , I XB (8 50 ) , I YB ( 350 ) 
*. I ZB ( 3 50 ) • I XY ( 850) ,IXZ(850) .IYZI850) ,CnSD(3,3) ,PI1 , P 1 2 , P 13 , SWE I GH, 
SXPAR.YRAR.ZBAR 

COMMON / T ON SOR/ A 1 1 . A 12 , A 13 . A 22 , A2 3, A33 
ECU IVAL ONCE ( A 12. A21 ) , ( A 13. A31) , ( A2 3 , A3 2) 

REAL (XY, IXZ. I YZ.MASS .! XC , IYC. IZC. 1X8,1 YB.I ZB 
WRI TF( 6. 2 13 ) 

2 1 3 FORM A T( 1 H 1 ) 

WRITE (6, 216) 

21 6 FORMA r< 1 HL ) 

WR I TE( A, 2 15 ) 

2(5 FORMAT ( ,30X. 46H *■*"*'*** COMPONENT MASS CHARACTER 1ST I CS ******) 

WR I TE! 6. 2 16) 

C 

C INDIVIDUAL COMPONENT OLTPUT BLOCK 

C 

N= ICNT 
OCNT-O 
or 6 K= 1 ,N 
J CNT =.) CNT + 1 

WR ITF (6.200) JCNT • ( Of .SCR I ( K , J ) , J= 1 . 5 ) , NCCM P ( K ) , I PART ( K ) , MO DEL ( K ), 
*ICATF:(K) 

200 FORMAT! 2X, DELIST NO. . I 5 ,4X .5A6.9HC0MP. N0..I6.10H PART NO. ,16, 
l L I H MODEL NO. , I6.L9H 0AT6 OF ENTRY- ,16) 

WRITE! 6, 2 01 ) MA SSI K ) , XC ( K ) ,YC ( K) ,ZC ( K) 

201 FORMAT! 6X. 8HWF IGHT= .E12.4.18H POSITION- XC=,E11.4, 

*9E INCHES .6X.6H YC=,EI1.4,9H INCHES ,6X,AH ZC=,E11.4,8H INCHF 

*S ) 

WRITE! 6, 251 ) I XC ( K 1 , 1 YC ( K ) , I ZC ( K ) 

251 FORMAT! 6X. 38HTNPUT MOMENT OF INERTIA- IXX=,E11.4, 

*9H SL UG- F T2 , 8X , AH I Y Y- ,E .1 l . 4 , 9H SLUG-F T2 ,6X , AH I ZZ= , El 1 . 4 , 9H SUJG-FT 

* 2 ) 

WRITF16.25?) PXYOIK) ,PXZC(K) ,PYZO(K) 

252 FORMAT! 6X, 38HINPUT PRODUCT OF INERTIA- PXYO=, El 1 . A, 

*9h SL UG- F T 2 , 7X , 5HP XZ 0-= , E 1 1 . A, 9H SLUG-F T2 ,5 X , 5H P YZO= , El 1 . 4 ) 

WRITE- (6, 214) XWR TCG! K) ,YWR ICG ( K) , ZWRTCG ( K) 

214 FORMAT! 38H COMPONENT POSITION WRT SYSTEM CM,7H CM-X= .E10.3, 

1 9 E INCHES , 6 X , 7H CM-Y= ,E10.3,9H INCHES ,6X,5HCMZ= ,E10.3, 

19 F INCHES ) 

WRITE! 6, 212) I XB ( K ) , I YB ( K ) , IZB !K) 

212 FORMAT (3BH MOMENT OF INERTIA CONTRI BUT I ON - , 7 H IXC= ,E10.3, 

1 9 E SLUG-F I 2,6X, 7H IYC= ,E1C.3,9H SLUG-F T2 ,6X ,5HI Z C= .E10.3, 

1 9 E SL UG- F T 2 ) 

WRITEI6.234) I X Y (K ) . I X Z ( K ) ,IYZ(K) 

234 FORMAT! 3811 PRODUCT OF INERTIA CONTR I BUT I ON-, 7H PXY= , E10.3, 

1 9 E SLUG— F T 2, 6X , 7H PXZ= ,E1C.3,9H SLUG-F T2 ,6X , 5HPY Z= ,E10.3, 



r> n non 


19F SLUG-FT2) 

WR I TE( 6, 2 1 1 ) 

21 l FORMAT ( 1HL ) 

6 CONTINUE 

INERTIA TENSOR OUTPUT BLOCK 

WRI TE! 6, 2 13 ) 

WRITE! 6. 2 17) 

2 1 7 FORMAT ! 1H4 ) 

WR ITE( 6, 2 18 ) 

218 FORMAT ( 40X.48HINERTIA TENSOR AT SYSTEM CENTER OF MASS 
WRITE! 6. 2 16) 

WRITE(6.219) 

219 FORMAT! 1 5X , 2H* * , 2 3X , 2H** . 1 1 X ,2 H** , 46 X . 2H** ) 

WR ITF! 6, 220 ) 

220 FORMAT! 15X.2H* .23X.2H *,UX,2H* .46X.2H *) 

WRITE! 6.221) All.A12.A13 

221 FORMAT ( .1 5 X « 6 H* I 1 1 . 6X , 3H 1 1 2 , 6 X , 6H I 1 3 *,llX,2h* .EU.4.6X, 

IE 1 1 .4. 6X.E1 1.4. 3H *) 

WRITE! 6.220) 

WR ITE(6.222 ) A21.A22.A23 

222 FORMAT! 15X. 6H* 1 2 1 . 6X ,3H I 22 ,6 X , 6H I 23 * ,5 X , 1 H= , 5 X , 2H* .Eli. 4, 

16X, El 1 .4.6X.E1 1.4.3H *) 

WRITE! 6. ? 20 ) 

WRITE! 6.223) A21.A32.A33 

223 FORMAT! 15X. 6H* I 3 1 . 6X , 3H I 32 ,6 X ,6H I 33 *,11X,2H* .FI1.4.6X, 

1E1.1 .4.6X.E 1 1.4, TH *) 

WRITE! 6, 220) 

WRITE! 6.219) 

WRITE! 6. 2 17 ) 

WRITE! 6, 224 ) 

224 FORMAT! 40X. 37HPRINCIPAL MASS MOMENTS OF INERTIA) 

WR ITF! 6, 2 1 1 ) 

WRITE! 6 « 232 )PI I.PI2.PI3 

232 FORMAT ( 8H0R00 Tl= .E12.5.9H SLUG-E T2 ,6X , 7UROOT2= .E12.5-.9H SLUG-FT2 
l, 6X. 7HRCI0T3 = .E12.5.9H SLUG-FT2) 

DIRECTION COSINE OUTPUT BLOCK 

WR ITF! 6,225) 

225 FORMAT! 1 HL ) 

WRITE! 6, 226) 

226 FORMAT! 40X. 32MOIRECTION COSINE ROTATION MATRIX) 

WRITE! 6,216) 

WR ITF! 6, 227 ) 

WRITE(A,228) 

227 FORMAT! 25X , 2H* * . 21 X , 2H** , l 3X.2H** .26X.2H**) 

228 FORMAT! 25X, 2H* .21X.2H *,13X,2H* .26X.2H *) 

WRITE! 6,229) (COSL) ! 1 . J ) . J = 1 ,3) 

229 FORMAT (25X.7H* L 1X8 , 4X , 4HL 1 YB ,4 X ,bHL 1 /.B *,13X,3H* .F7.4.1X, 

1E7.4, IX, F7.4, 3X, 1H* ) 

WRITE(6,223) 

WRITE! 6,230) (C0SD12.J) .J=l,3) 

2 30 FORMAT! 25X. 7H* L2 XB.4X.4HL2YB .4X.6HL2ZB * ,6X , 1H= , 6X ,3H* 

1F7.4, IX.F7.4.LX.F7.4.3X, Lit*) 
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WRITF(A.228) 

WRITE) 6. 231 ) CCnSO ! 3, J.) ,J = 1 ,3) 

23 1 FORM AT ( 2 5X • 7H* L3 XB , AX , 4HL 3YB , AX . 6HL3 ZB *,13X,3H* .F7.A.1X, 

IF 7. A. 1X.F7.A, 3X » 1H* ) 

WRITE<6.228) 

WRITE! 6. 227) 

C 

C SYSTEM OUTPUT BLOCK 

r. 

W R I T F ( 6 . 2 1 1 ) 

WR ITE! 6. 209) SWFIGH 

20<5 FORMAT ( 16H SYSTEM WEIGHT= .E12.5.7H POUNDS) 

WR ITE! 6 . 2 .3 6 I 
WRITE! A, 205) 

20 5 FORMAT ! 30X , A5HSYSTEM CENTER OF MASS WR T REFERENCE SYSTEM * 
WRITF! 6.206) XB AR , Y8 AR . ZBAR 

206 FORMAT! 7E0XBAR= .E12.5.7H INCHES. 6X.7H YBAR= .E12.5.7H INCHES, 
16X.7H ZBAR = .F12.5.7H INCHES) 

WR I TE< 6, 236 ) 

W R I T F ( 6 , 2 0 8 ) N 

208 FORMAT ( 29H TOTAL NUMBER OF CCM'PCNENTS= ,16) 

236 FORMAT! LHJ ) 

R FTtJRN 
END 
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SIRFTC F I GEN 


SIRROUTINE E I GEN < EVL .EV2 ,EV3) 

EGUIVALENCE(A12.A21 ) , ( A31 ,A13) , IA23.A32) 

GCMMON/ TENSOR/ A1 1, A 12, A 13 ,A22.A23,A33 
CEMMON/CUEFF/P «0«R 
P=-( A1 1+A22+A33) 

0=-(-A11*A2 2-A11*a3 3— A22*A33+A 12*A21+A13* A 31 +A32*A23 ) 

R ='— ( AL A 22* A 3 3- A 12* A 21* A 3 3 + A l 2*A31*A23+ A1 3*A2 l * A3 2- A1 3* A2 2*A31 
1- A3 2 * A 2 3 * A 1 1 ) 

CALL CUBIC ( EVl .EV2»EV3) 

RFTURN 

END 
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$ I ft FTC CUBIC 


S USROt : I INC CUBICIKI.R2.R3) 

CCMMON/COEFF/P .Q.R 
DATA CONRTD/57. 2957795/ 

A=( 1 .0/3.0 }*< 3. 0*G-P*P ) 

R = ( 1.0/ 27.0)* ( 2.0* ( P** 3 ) -9. 0*P * 3*2 7 • 0*R ) 

I F ( ( B* R / 4 . 0 + A* A * A / 2 7 . U ) .GT.O.O) GO TO 9 
I F( ( B*B/4.0+A*A*A/27.O).EO.0.O) GO TO 7 
C TFREF REAL AND UNEQUAL ROOTS 

P E I = AR CO S ( — B / ( 2.0* SORT <-< A* A* A ) /27.0) ) ) 

FAC TOR = 2 . 0*S0R T(— A/3.0) 

R 1 = F AC T UR*CO S ( PHI/3.0)— P/3.0 

R 2 = F A C T 0 R * C 0 S ( < PHI / 3. 0 ) + I 1 20. 0 /C ONRTD )) -P/3 . 0 
R 3=FACT0R*C0S( ( PHI / 2 . 0 ) * < 24 0. 0 /C CNR TD > >-P/3.0 
GC TO 5 

C TFRFE REAL ROOTS. AT LEAST TWO EQUAL 

7 AA=+(-B/2 .0 )**( 1 .0/3.0) 

R 1 = 2. 0*A A- P /3 . 0 
R ?=-AA -P/3.0 
R 3 =— A A— P / 3 . 0 
GC TO 5 

C ONE REAL ROOT , AND TWO CONJUGATE IMAGINARY 

9 R 1=0.0 
R 2=0.0 
R 3=0.0 
5 RETURN 
END 



APPENDIX D 


SAMPLE CALCULATION 

A sample problem and its solution are presented to demonstrate the preparation of 
input data and to illustrate the output format. Using the proper matrix relations, a check 
is made on the computed principal mass moments of inertia and the directions of the 
principal axes. 


Problem 


A composite body consists of six regular geometric shaped components located 
within a cylindrical shell. A top view of the body is presented in figure 5. Assume a 
reference system coincident with the geometric axes through the base of the shell. Com- 
ponent II has its principal axes rotated 45° from the reference axes. The principal mo- 
ments of inertia of component II are 

I x , = 2. 5 (slug-ft 2 ) 

I yt = 2. 5 (slug-ft 2 ) 

I z , = 1.0 (slug-ft 2 ) 

The inertia tensor of component II must be found about axes parallel to the reference 
axis system. The relation between the principal moments of inertia and moments of in- 
ertia about axes parallel to the reference system is 

[i] = [*] T [r] [;] (Di) 


where [Z ] is the direction cosine matrix between the principal axes and axes parallel to 
the reference system. For this case, [Z ] represents a single rotation in the xy plane 
through a, or 


m 


cos a 
sin a 


-sin a 
cos a . 


(D2) 


and 
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(D4) 


Substituting equations (D2), to (D4) into equation (Dl), then 

n r 2 2 

I xy I x , COS a +1 , sin a (I , - I x ,)sin a cos a 

(D5) 

2 2 

Iyx Iy (Iy, - I x t)sin a cos a 1^ sin a + 1 , cos a 
For component II, a = 45°, so equation (D5) becomes 



The input data required for the problem is shown in table I. Table n shows the input 
card requirements. 

The computer output for the problem is given in table III. The component mass 
characteristics are individually listed. The component moment of inertia contribution 
and product of inertia contribution represent the inertia of the component about the cen- 
ter of mass axis system. The inertia tensor is consistent with the definition of equa- 
tion (6). The principal mass moments of inertia are numbered to correspond with the 
nomenclature of the direction cosine matrix defined by equation (13). In the output direc- 
tion cosine matrix, L1XB represents the direction cosine between the axis about which 
ROOT1 is defined and the x axis of the center of mass coordinate system. The center 
of mass location is with respect to the reference axis system. 

The matrix [Z] in equation (Dl) is an orthogonal transformation and therefore pre- 

T 

multiplying (Dl) by [Z] and postmultiplying by [Z] 

[rj = mmm T m 


From matrix theory, the trace of any square matrix is invariant under an orthogonal 
transformation. By inspection of the computer results, the trace of [I] equals the trace 


27 



of [I 1 ]. As a further verification of the computer results, the matrices are substituted 
into equation (D6). 


0.2015 0.7473 

. 8440 . 1956 

. .4970 -.6349 


-0. 6332' 
.4994 


295. 8 


-5. 472 


-. 5915J 


L- 25. 32 


-5.472 

-25. 23' 


'0. 2015 

0. 8440 

0. 4970' 

343.0 

-15. 61 


.7473 

. 1956 

-. 6394 

-15. 61 

328. 3 . 


.-. 6332 

.4994 

-. 5915. 


'354. 754 
.013 
1.083 


0.013 
279. 585 
-.014 


- 1 . 11 ' 
-. 252 
333. 704_ 


TABLE I. - SAMPLE PROBLEM INPUT DATA 


Component 

Description 

Weight, 

lb 

Center of 
mass, 
ft 

2 

Mass moment of inertia, slug-ft 

s 



I 

Cylindrical shell: 
Length, 10 ft 
Inner radius, 24 in. 
Outer radius, 25 in. 

640 

(0,0, 5) 

208.4 

208.4 

83.4 

n 

Cube: 

1 by 1 by 2 ft 

160 

(1, b7) 

1.3 

1. 3 

1 

m 

Sphere: 

Radius, 1 ft 

128 

(hi 3 ) 

1. 6 

1. 6 

1. 6 

IV 

Thin plate: 
1 by 1 ft 

8 

('!’ ~2’ 4) 

0. 1 

0. 1 

0. 1 

V 

Cylinder: 

Length, 5 ft 
Radius, 1/2 ft 

80 


5.4 

5.4 

0. 3 

VI 

Cylinder: 

Length, 3 ft 
Radius, 1/2 ft 

64 

(-1, 1, 5) 

1. 6 

1.6 

0. 3 

vn 

Cylinder: 

Length, 8 ft 
Radius, 1/2 ft 

1024 

(-i i°’ 4 ) 

4 

' 

4 

174. 6 
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TABLE m. - SAMPLE PROBLEM OUTPUT LISTING 


****** COMPONENT MASS CHAR ACT ERl ST ICS ****** 


LIST NO. 1 CYLINDRICAL SHELL 

COMP. 

NO. 1 

PART NO. 

1 MODEL 

NO. 1 

OATE l 

OF ENTRY- 

- 

90268 

WEIGHT= 0.6400E 03 POSITION- XC= 

0 . 


INCHES 

YC = 

0 . 


INCHES 

7.C = 

0.6000E 

02 

INCHES 

INPUT FOMENT OF INERTIA- I XX= 

0.2084E 

03 

SLUG-FT2 

I YY = 

0.2084E 

03 

SLUG-FT2 

III- 

0. 8340E 

02 

SLUG-FT2 

INPUT PRODUCT OF INERTIA- PXYO= 

0 . 


SLUG-FT 2 

PX 2 0= 

0 . 


SLUG-FT2 

PYZO = 

0 . 



COMPONENT POSITION WR T SYSTEM CM CM-X” 

-0.741E 

01 

I NCHES 

CM— Y = 

0.1 16E 

01 

INCHES 

CMZ” 

-0.659E 

01 

INCHES 

MOMENT OF INERTIA CONTRIBUTION - I XC= 

0.215E 

03 

SLUG-FT2 

* IYC= 

0.222E 

03 

SLUG-FT2 

I zc = 

0.912E 

02 

SLUG-FT2 

PRODUCT OF INERTIA CONTRIBUTION- PXY= 

-0.1 19E 

01 

SLUG- FT 2 

PXZ = 

0.675E 

01 

SL UG-FT2 

PYZ = 

-0. 106E 

01 

SLUG-FT2 


LIST NO. • 2 CUBE 


COMP. 

NC. 2 PART 

NO. 

1 MOOEL 

NO. 1 

DATE i 

OF ENTRY- 


9026 8 

WEIGHT* 0. 1600E 03 POSITION- XC= 

0. 

1200E 

02 

INCHES 

YC = 

0.1200E 

02 

INCHES 

ZC- 

0.8400E 

02 

INCHES 

INPUT MOMENT OF INERTIA- I XX= 

0. 

1300E 

01 

SLUG-F T2 

I YY= 

0.I300E 

01 

SLUG-FT2 

I ZZ* 

0. 1000E 

01 

SLUG-FT2 

INPUT PROOUCT OF INERTIA- PXYO= 

0. 



SLUG-FT2 

PXZO= 

0. 


SLUG-FT2 

pyzo = 

0. 



COMPONENT POSITION WR T SYSTEM CM CM-X* 

-o 

.1946 

02 

I NCHES 

CM — Y = 

-0. 108E 

02 

INCHES 

CMZ ” 

-0.306E 

02 

INCHES 

MOMENT OF INERTIA CONTRIBUTION - ! XC= 

0 

• 3 77E 

02 

SLUG-FT 2 

IYC = 

0.4666 

02 

SL UG- F T2 

I ZC- 

O.lBtE 

02 

SLUG-FT2 

PROOUCT OF INERTIA CONTRIBUTION- PXY= 

0 

.727E 

01 

SLUG-FT2 

PXZ = 

0.205E 

02 

SLUG-FT2 

PYZs 

0. 114E 

02 

SLUG-FT2 


LIST NO. 3 SPHERE 

CCMF. 

NO. 3 

PART NO. 

I MODEL 

NO. 1 

DATE OF ENTRY- 


90268 

WFIGHT- 0.1280E 03 POSITION- XC” 

0.6000E 

01 

l NCHES 

YC= 

0.60 OOE 

01 

INCHES 

ZC= 0. 3600E 

02 

INCHES 

INPUT MOMENT OF INERTIA- I XX= 

0. 16006 

01 

S LUG-FT2 

IYY= 

0.1600E 

01 

SLUG-F T2 

I 17= 0. 1600E 

01 

SLUG-FT2 

INPUT PRODUCT OF INERTIA- PXYO= 

0 . 


SLUG-FT 2 

PXZO= 

0 . 


SL UG-FT2 

PYZO = 0. 



COMPONENT POSITION WKT SYSTEM CM CM-X= 

-0.134E 

02 

INCHES 

CM — Y = 

-0.484E 

01 

INCHES 

CMZ= 0.1 74E 

02 

INCHES 

MOMENT OF INERTIA CONTRIBUTION - I XC= 

0.1 06E 

02 

SLUG-FT 2 

I YC= 

0.149E 

02 

SLUG-FT2 

IZC= 0.722E 

01 

SLUG-FT2 

PRODUCT OF INERTIA CONTRIBUTION- PXY= 

0.179E 

01 

SLUG-FT 2 

PXZ = 

-0.645E 

01 

SLUG-FT2 

PYZ * -0.233E 

01 

SLUG-FT2 


LIST NO. 4 THIN PLATE 

CGMP. 

NO. 4 

PART NO. 

l MODEL 

NO. 1 

OATE OF ENTRY 

- 

90268 

HEIGHT- 0.8000E 01 POSITION- XC*- 

-0.6000E 

01 

INCHES 

YC=- 

0.6000E 

01 

INCHES 

ZC= 0.66006 

02 

INCHES 

INPIIT MOMENT OF INERTIA- I XX= 

0. 10006 

00 

SLUG-F T2 

I YY= 

0 . 10 OOE 

00 

SLUG-F T2 

IZZ= 0.1000E 

00 

SLUG-F T2 

INPUT PRODUCT OF INERTIA- PXYO” 

0. 


S LUG-FT2 

PXZO= 

0. 


SLUG-FT2 

P YZO = 0. 



COMPONENT POSITION WR T SYSTEM CM CM-X= 

-0.141E 

01 

INCHES 

CM — Y- 

0.716E 

01 

INCHES 

CMZ = -0.126E 

02 

INCHES 

MOMENT OF INERTIA CONTRIBUTION - I XC= 

0.4626 

00 

SLUG-FT2 

I YC= 

0.377E 

00 

SLUG-F T2 

IZC= 0.192E 

00 

SLUG-FT 2 

PRODUCT OF INERTIA CONTRIBUTION- PXY= 

—0. 1 75E-01 

SLUG-FT2 

PXZ = 

0.308E- 

01 

SL UG- F T2 

PYZ= -0. 156E 

00 

SLUG-FT2 


LIST NO. 5 CYLINDER L=5.R=l/2 


COMP. 

NO. 5 

PART NO. 

1 MODEL 

NO. 1 

OATE l 

OF ENTRY- 


90268 

WE IGHT = C.8000E 02 POSITION- 

XC = 

0. 12 OOE 

02 

I NCHES 

YC=- 

0.1200E 

02 

INCHES 

ZC = 

0.3000F 

02 

INCHES 

INPUT MOMENT OF INERTIA- 

1 xx = 

0.5400E 

01 

S LUG-FT2 

I YY- 

0.5400E 

01 

SLUG-F T2 

!ZZ = 

0.3000E- 

00 

SLUG-FT2 

INPUT PRODUCT OF INERTIA- 

PX Y0= 

o. 


SLUG-FT2 

PXZO = 

0. 


SLUG-FT2 

PYZO = 

0. 



COMPONENT POSITION WRT SYSTEM CM 

CM-X= 

-0.194E 

02 

INCHES 

CM ~Y= 

0.132E 

02 

INCHES 

CMZ = 

0. 2 34E 

02 

INCHES 

MOMENT OF INERTIA CONTRIBUTION - 

I XC = 

0.1 79E 

02 

SLUG-FT2 

IYC= 

0.214E 

02 

SLUG-FT2 

IZC* 

0.980E 

01 

SLUG-FT2 

PRODUCT OF INERTIA CONTRIBUTION- 

PXY = 

-0.441F 

01 

SLUG-FT2 

PXZ = 

-0.785E 

01 

SLUG-F T2 

PYZ = 

0.532E 

01 

SLUG-FT? 


LIST NO. 6 CYLINDER L = 3,R = l/2 

COMP. 

NO. 6 

PART NO. 

1 MOOEL 

NO. 1 

DATE OF ENTRY- 


90268 

WEIGHT* 0.6400E 02 POSITION- XC= 

-0.1200E 

02 

INCHES 

YC= 

0.1200 E_ 

02 

INCHES 

ZC= 0.6Q00E 

02 

INCHES 

INPUT .nUMENT OF INERTIA- I XX= 

O'. 1600E 

01 

SLUG-FT2 

[ YY= 

0. 16 00E~ 

01 

SLUG-F T2 

IZZ= 0.3000E 

00 

SLUG-FT2 

INPUT PRODUCT OF INERTIA- PXYO= 

0 . 


SLUG-FT 2 

PXZO* 

0 . 


SLUG-FT2 

P YZO= 0. 



COMPONENT 'POSITION WRT SYSTEM CM CM-X= 

0.459E 

01 

INCHES 

CM -Y= 

-0.108E 

02 

INCHES 

CMZ* -0.659E 

01 

INCHES 

MOMENT OF INERTIA CONTRIBUTION - I XC= 

0.382 E 

01 

SLUG-FT 2 

I YC= 

0.249E 

01 

SLUG-F T2 

IZC= 0.221E 

01 

SLUG-FT2 

PRODUCT OF INERTIA CONTRIBUTION- PXY= 

-0.686E 

00 

SLUG-FT2 

PXZ = 

— 0.4 18 E 

00 

SLUG-FT2 

P YZ = 0.987E 

00 

SLUG-FT2 


LIST NO. 7 CYLINDER L=8,R=l/2 

CCMF. 

NC. 7 

PART NO . 

1 MODEL 

NO. 1 

OATE OF ENTRY- 


90268 

WEIGHT” 0.1024E 04 POSITION- 

XC=-0. I300E 

02 

INCHES 

YC= 

0 . 


INCHES 

ZC= 0. 48 OOE 

02 

INCHES 

INPUT MOMENT OF INERTIA- 

1 XX= 0.4000E 

01 

SLUG-F T2 

I YY= 

0 .40 00 E 

01 

SLUG- FT2 

IZZ* 0. 1 746 E 

03 

SLUG-FT2 

INPUT PRODUCT OF INERTIA- 

PX Y0= 0. 


SLUG-FT2 

PXZO= 

0. 


SLUG-FT2 

PYZO= 0. 



COMPONENT POSITION WRT SYSTEM CM 

CM-X= 0.L06E 

02 

INCHES 

CM — Y = 

0.116E 

01 

INCHES 

CMZ = 0.541E 

01 

INCHES 

MOMENT. OF INERTIA CONTRIBUTION - 

I XC= 0.1C8E 

02 

SLUG-FT 2 

I YC= 

0.352E 

02 

SLUG-FT2 

1ZC= 0.200E 

03 

SLUG-F T2 

PRODUCT OF INERTIA CONTRIBUTION- 

P XY= 0.272E 

01 

S LUG- FT 2 

PXZ = 

0.126E 

02 

SLIJG-FT2 

P YZ = 0.139E 

01 

SLUG-FT2 


L 1ST NC. 8 END OF DATA 

CC-MP. NC. 0 

PART NO. 

0 MODEL 

NO. 0 

DATE OF ENTRY- 

0 

WEIGHT” 0. POSITION- XC= 

0. INCHES 

YC = 

0. 

INCHES 

ZC= 0. 

INCHES 

INPUT MOMENT OF INERTIA- I XX* 

0. SLUG-FT 2 

I YY= 

0. 

SL UG-FT2 

IZZ = 0. 

SLUG-FT2 

INPUT PRODUCT OF INERTIA- PXYO= 

0. SLUG-FT2 

PX Z0= 

0. 

SL UG-FT2 

PYZO* 0. 


COMPONENT POSITION WRT SYSTEM CM CM-X= 

-0.741E 01 INCHES 

CM— Y= 

0.I16E 01 

INCHES 

CMZ* 0.534E 02 

INCHES 

MOMENT OF INERTIA CONTRIBUTION - I XC= 

0. SLUG-FT 2 

IYC= 

0. 

SLUG-FT2 

IZC = 0. 

SLUG- FT2 

PRODUCT OF INERTIA CONTRIBUTION- PXY= 

-0. SLUG-FT2 

PXZ = 

-0. 

SLUG-FT2 

PYZ” 0. 

SLUG-FT2 


30 



TABLE B3. - Concluded. SAMPLE PROBLEM OUTPUT LISTING 

INERTIA TENSOR AT SYSTEM CENTER OF MASS 


** 




*<■ 


*• 






** 

* 




p 


* 






* 

* 

111 

112 

1 13 

* 


* 0.2958E 

03 

-0.5472E 

01 

— 0.2523E 

02 

* 

p 




* 


* 






* 

* 

I?t 

122 

I 23 

* 

= 

* -0.5472E 

01 

0.3430E 

03 

-0.1561E 

02 


* 




* 


* 






* 


131 

1 3? 

133 

* 


* -0.2523E 

02 

-0.1561E 

02 

0.3283E 

03 

* 

* 




* 


* 






* 

pp 




p p 


*# 






** 


PRINCIPAL MASS MOMENTS OF INERTIA 


ROOT 1= 0.35479E C 3 SLUG-F T2 RCOT2 = U.27959E 03 SLUG-FT2 R00T3= 0.33279E 03 SLUG-FT2 

DIRECTION COSINE ROTATION MATRIX 





** 


pp 




pp 

* 



* 


p 




p 

* 

1 1XR 

L1YB 

L1ZR * 


p 

0.2015 

0.7473 

-0.6332 

p 

* 



* 


p 




P 

* 

L2XR 

L2YB 

L2ZB * 

- 

p 

0.8440 

0. 1956 

0.4994 

P 

* 



p 


p 




p 

* 

L3XR 

L3YB 

L3Z0 * 


p 

0.4970 

-0.6349 

-0. 5915 

P 

* 



p 


p 




P 

** 



pp 


pp 




pp 


SYS T FM WEIGHT= C.21040E 04 POUNDS 

SYSTEM CENTER OF MASS WRT REFERENCE SYSTEM 

XBAR = -C.74144E 01 INCHES YRAK= 0.11635E 01 INCHES ~lb AK= 0.53407E 0^ INCHES 

TOTAL NlJVHER OF CCMPCNENTS= 8 
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Figure 3. - INRTIA calculates principal moments of inertia and direction of prin- 
cipal axes. 
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Figure 4. - OUT prints mass characteristics 
of each component and entire system. 
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